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ABSTRACT 

There is much evidence to suggest that stellar wind capture, rather than 
Roche lobe overflow, serves as the accretion mechanism onto the compact sec- 
ondary object in the massive X-ray binary LS 5039. The lack of significant 
emission combined with only a modest X-ray flux provide observational evidence 
that no large-scale mass transfer is occurring (consistent with our estimate of 
the radius of the 06.5 V((f)) optical star that is smaller than its critical Roche 
radius). Here we determine the mass loss rate of the optical star from the broad, 
residual emission in the Ha profile. Using a stellar wind accretion model for a 
range in assumed primary mass, we compute the predicted X-ray luminosity for 
the system. We compare our results to the observed X-ray luminosity to deter- 
mine the mass of the compact object for each case. The companion appears to 
be a neutron star with a mass between 1 and 3Mq. With our new constraints on 
the masses of both components, we discuss their implications on the evolution 
of the system before and after the supernova event that created the compact 
companion. The binary experienced significant mass loss during the supernova, 
and we find that the predictions for the resulting runaway velocity agree well 
with the observed peculiar space velocity. LS 5039 may be the fastest runaway 
object among known massive X-ray binaries. 

Subject headings: binaries: spectroscopic — stars: early-type — stars: individual 
(LS 5039, RX J1826.2-1450) — X-rays: binaries 
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1. Introduction 

LS 5039 is a relatively faint {V = 11.3) and massive star of type 06.5 V((f)) (Clark et al. 
2001) that is one of only a few confirmed massive X-ray binaries (MXRBs) with associated 
radio emission (Ribo et al. 1999). It has radio-emitting relativistic jets characteristic of 
galactic microquasars, and it is probably a high energy gamma ray source as well (Parades 
et al. 2000). We recently discovered that the system is a short period binary (P = 4.117 ± 
0.011 d) with the highest known eccentricity (e = 0.41 ± 0.05) among O star binaries with 
comparable periods (McSwain et al. 2001a). This high eccentricity probably results from 
the huge mass loss that occurred with the supernova (SN) explosion that gave birth to the 
compact star in the system (Bhattacharya & van den Heuvel 1991; Nelemans et al. 1999). 
Binaries that suffer large mass loss in a SN are expected to become runaway stars, and 
recently both we (McSwain et al. 2001b) and Ribo et al. (2002) found that LS 5039 has a 
relatively large proper motion that indicates that the binary has a record-breaking peculiar 
space velocity among MXRBs. 

Reliable estimates for the masses of the components of the binary are of key importance 
for any discussion about the evolution of this remarkable system. Here we present an inves- 
tigation of the possible mass range for the X-ray star based on a stellar wind accretion model 
for the X-ray production (§2). This analysis relies on our observations of the wind emission 
effects in the Ho; profile and the wind models of Puis et al. (1996). We then apply our derived 
mass estimates for both stars to determine the probable masses and orbital parameters prior 
to the SN based on the system's current eccentricity (§3). We find that the predicted and 
observed runaway velocities are in good agreement, and, thus, LS 5039 provides the best 
verification to date of model predictions about the outcome of a SN explosion in a massive 
binary. 

2. X-ray Fueling by Wind Accretion 

The X-ray production in MXRBs results from mass accretion through a Roche lobe 
overflow (RLOF) stream or by Bondi-Hoyle accretion from the stellar wind of the luminous 
primary (Kaper 1998). The systems experiencing RLOF usually have large X-ray luminosi- 
ties and striking optical emission lines. However, LS 5039 has a modest X-ray flux (Ribo et 
al. 1999) and no obvious emission lines (Clark et al. 2001; McSwain et al. 2001a), and we 
show below that the 0-star is probably smaller than its critical Roche surface. Thus, we 
can safely assume that most of the mass accretion in LS 5039 occurs through capture of the 
primary's stellar wind flow. 
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The wind accretion luminosity of MXRBs depends on the system separation, wind 
velocity law, mass loss rate, and mass of the accretor (Lamers et al. 1976). If we can estimate 
these binary and wind parameters from spectroscopy, then the mass of the compact star can 
be found by comparing the predicted and observed X-ray luminosities, Lx- Here we present 
such an analysis based on the mass loss rate derived from our observations of the Ha profile 
(McSwain et al. 2001a). All the parameters needed in this analysis ultimately depend on the 
assumed mass of the 0-star primary that is poorly constrained at the moment. The simplest 
assumption is that the primary has a mass typical of single stars of its spectral classification, 
approximately AQMq (see the mass calibration of Howarth & Prinja (1989) and the study of 
comparable stars in the young binary, DH Cep, by Penny et al. (1997)). On the other hand, 
there is evidence that the primaries in MXRBs may be undermassive for their luminosity 
(Kaper 2001), and in extreme cases, their mass may be a factor of three lower than the mass 
derived by comparing their position in the Hertzsprung-Russell diagram with evolutionary 
tracks. Thus, we show here the results of a wind accretion model for LS 5039 based on a 
range in assumed primary mass of 10 - 40 Mq. 

Puis et al. (1996) show how the Ha profile in 0-stars grows from a pure absorption 
feature to a strong emission line with increasing stellar wind mass loss rate, and they present 
a scheme to estimate the mass loss rate based on the observed equivalent width, of Ha. 
We observed the Ha hne in the spectrum of LS 5039 over three runs in 1998, 1999, and 
2000 with the Coude Feed Telescope at Kitt Peak National Observatory (see McSwain et 
al. (2001a) for details), and we show in Figure 1 the average profile from each run after 
shifting each spectrum to the rest frame and convolving the 1998 and 2000 spectra to the 
resolution of the 1999 run (i? = A/AA = 4000). The Ha profile appears to be filled in by 
broad, residual emission that appears as shallow emission peaks in the wings (perhaps also 
present in the vicinity of He I A6678 and He II A6683). The emission was apparently stronger 
during the 1998 and 1999 runs, and we compare in Figure 1 these profiles with the deeper 
absorption observed in 2000 {dotted lines). The net equivalent width over the entire emission 
and absorption blend was 2.22, 2.15, and 3.10 A (±10%) for the average profiles from 1998, 
1999, and 2000, respectively. We use the mean of these values, 2.49 A, in the analysis below. 
We see no evidence in our spectra of variations in the residual emission with orbital phase, 
and so we can reliably assume that this weak emission forms in the wind of the primary 
(and not, for example, in a gas stream or accretion disk in which we would observe orbital 
variations in radial velocity). 

We need estimates of the underlying photospheric component of absorption equivalent 
width, effective temperature, Tefj, radius, Rq, and wind velocity law in order to derive a mass 
loss rate from the observed equivalent widths using the wind models of Puis et al. (1996) 
(see their Table 7). We adopted the spectral classification calibration of Howarth & Prinja 
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(1989) to estimate T^s = 40.5 ± 2 kK and log^ = 4.0 ± 0.2 for the 06.5 V((f)) primary, 
which is in reasonable agreement with the results from Stromgren photometry by Kilkenny 
(1993) and with the values of these parameters derived by Puis et al. (1996) for similar 
stars. The radius, Rq, follows from logg and the assumed mass, Mq- We assumed the wind 
velocity law usually applied to 0-stars: v{r) = foo(l — Roh'Y with {3 = 0.8, foo = 2.6fcsc, 
Vesc = (2(1 - ^)G'Mo/i^o)^/^ and T = 2.6 x IQ-' {Lq / Lq)/ (Mq/Mq) (Howarth & Prinja 
1989; Puis et al. 1996; Lamers & Cassinelli 1999). The final parameter to be set is the 
appropriate photospheric absorption equivalent width. Unfortunately, Puis et al. (1996) only 
give one value of Wx — 3.29 A that is specified for a model photosphere with Tes — 40 kK 
and loggi = 3.7. While this effective temperature is nearly applicable here, we expect the 
equivalent width will be larger in our logg = 4.0 case since the Stark broadening of the 
Balmer lines increases with gravity. We used the NLTE models of Auer & Mihalas (1972) 
to estimate that the Ha equivalent width should be 10% stronger at logg = 4.0 compared 
to logg = 3.7 (for Tes = 40 kK), and we pro-rated accordingly the photospheric equivalent 
width given by Puis et al. (1996). Our derived estimates of Rq, f oo, and mass loss rate log M 
are given for several assumed values of Mq in Table 1. The mass loss rate values are close 
to the average found for 06.5 V stars by Howarth & Prinja (1989), logM = —6.3 (in units 
of Moy-i). 

We next calculated the predicted X-ray luminosity for a range in assumed secondary 
mass, Mx, using the wind accretion model of Lamers et al. (1976) (using the wind velocity 
law noted above). Their expression for the X-ray accretion luminosity is 

Lx = {3x lO^-Uo) (S^Mx/Mq 

where the efficiency factor for the conversion of accreted matter to X-ray fiux is ^ ~ 0.1. 
The accretion rate is 

Sa = 7TrlM/{A7ra^) 
for a separation a and an accretion radius 

r„ = 2GMx/vl, 

where Vrei is the fiow velocity of the wind relative to the accreting star (see their eq. 10a). 
Given the two assumed masses and the known period, we determined the time-averaged 
separation of the binary and then calculated Lx using our estimates of the wind speed and 
mass loss rate from above. A sample result is shown in Figure 2 for the average and extreme 
values of the derived mass loss rate (from the 1999 and 2000 Ho; equivalent widths) for a 
test primary mass of Mq — 25Mq. 

The observed value of Lx depends on the assumed distance. We used the magnitude V — 
11.35±0.1 and reddening E{B-V) = 1.2±0.1 from Clark et al. (2001), and then estimated 
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the star's angular diameter by comparing the unreddened magnitude with model magnitudes 
from Kurucz (1994) for the adopted T^g and log (7 (calibrated with the fundamental Teg data 
of Code et al. (1976)). The distance is then found from the working value of the stellar radius, 
d = 0.32 Rq/Rq kpc. The observed X-ray luminosity was adjusted from the measurements 
of Ribo et al. (1999) who assumed a distance of 3.1 kpc. This prorated estimate of Lx is 
shown as the shaded region in Figure 2. The estimated mass of the secondary is found where 
the predicted wind accretion luminosity matches the observed Lx- 

We repeated the wind accretion calculation for a grid of assumed 0-star masses, and 
our final results arc listed in Table 1 and plotted in a mass diagram in Figure 3. The 
results are relatively insensitive to our assumptions about T^s and logs'. We found, for 
example, that models made using log^f — 3.7 led to lower mass loss rates but also lower wind 
velocities, so that the wind accretion rates changed very little. Because the slope dLx / dMx 
is relatively large (see Fig. 2) , the details of the accretion model (for example, the value of 
the efficiency parameter, Q do not greatly affect the implied secondary masses. In fact, the 
largest uncertainty in the results comes from the variation in Ha equivalent width between 
observing runs, which presumably refiects significant changes in the stellar wind mass loss 
rate. Table 1 also lists the size of the Roche radius at periastron, and we find that the O- 
star fits comfortably within its Roche lobe over the full range in assumed mass. Our results 
indicate that the secondary has a mass between 1 and 3M0 (based on the extremes of the 
observed Ha equivalent width and adopted primary mass range), and, thus, the secondary 
is probably a neutron star. 



3. A Supernova in a Bineiry 

The catastrophe of a SN explosion in a massive binary has two immediate consequences. 
First, the system acquires a large eccentricity that is directly related to the amount of 
mass lost in the SN event. The periastron separation in the altered orbit corresponds to 
the pre-supernova semi-major axis. Second, by conservation of momentum, we expect the 
entire system to attain a runaway velocity that again depends on the mass lost in the SN. 
We demonstrate here that the values of the observed eccentricity and the probable masses 
derived above indicate that a huge amount of mass was lost in the SN event that formed 
LS 5039. 

We can use the observed eccentricity to relate the pre- and post-SN orbital parameters 
if we make the following reasonable assumptions: (1) the pre-SN orbit was circular (almost 
certainly the case for such a short period and evolved system), (2) any kick velocity imparted 
to the remnant core due to asymmetries in the explosion was relatively small (Nelemans et al. 
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1999), (3) the primary suffered only minor ablation of mass in the explosion so that its pre- 
and post-SN mass is the same (Fryxell & Arnett 1981), and (4) the system has experienced 
little or no tidal reduction of the orbital eccentricity since the SN event. Bhattacharya & 
van den Heuvel (1991) and Nelemans et al. (1999) give the expressions required to calculate 
the pre-SN parameters, and wc list in Table 2 the resulting orbital and physical parameters 
for an eccentricity e = 0.41 ± 0.05 (McSwain et al. 2001a) and the masses derived in §2. 
Table 2 gives the pre-SN period, p^^^*^^^^ semi-major axis, a™'*''^^ SN precursor mass, Mg"'*'^', 
and the mass lost in the SN, AM2. 

The first striking result is the large mass loss that occurred in the SN in LS 5039. Some 
5 - 17 Mq was lost in the explosion, amounting to more than 81% of the precursor's mass. 
This is much larger than the typical inferred mass loss fraction of 35% found by Nelemans 
et al. (1999) for black hole X-ray binaries. The SN precursor had a smaller mass than the 
primary at the time of the explosion, so the system was not in danger of total disruption 
(Nelemans et al. 1999). 

The pre-SN orbit was very compact and the stars were in close proximity. We list in 
Table 2 the sizes of the Roche radii in the pre-SN stage. These radii are quite restrictive, 
and we find that the primary overfills the Roche lobe in the higher mass solutions. However, 
the current radius Rq given in Table 2 may be larger than the radius at the time of the 
SN if the primary has evolved to a larger size since then. Nevertheless, the secondary's 
Roche radius is also quite small, and it is possible that the system was in a contact or 
over-contact configuration at the time of the SN. Some kind of close interaction must have 
occurred since by the time of the SN the mass ratio had reversed and the separation had been 
increasing (Wellstein et al. 2001). Given this evidence of a pre-SN interaction, the observed 
C deficiency of the primary (McSwain et al. 2001a) probably results from nuclear-processed 
gas transferred from the SN progenitor. The currently faster than synchronous rotation of 
the primary (McSwain et al. 2001a) may correspond to the synchronous rate in the shorter 
period, pre-SN configuration. 

Finally, we list in Table 2 the predicted system runaway velocity, Vgys, based on conser- 
vation of momentum (with errors propagated from the uncertainty in the observed eccentric- 
ity). The predictions suggest that the system should have a runaway velocity in excess of 100 
km s~^. Although the systemic velocity along the line of sight is unexceptional (McSwain 
et al. 2001a), the tangential velocity does appear to be large. The system has a proper 
motion of fi^ = 4.8 ± 2.4 mas y""*^ and fi^ = —12.2 ± 2.3 mas y^^ in the Tycho-2 catalogue 
(Hog et al. 2000), and Ribo et al. (2002) have recently used optical and radio astrometry 
to find an improved estimate of fj^a = 4.7 ± 1.1 mas y~^ and fj^s = —10.6 ± 1.0 mas y~^. 
We used the latter measurement together with the systemic radial velocity from McSwain 
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et al. (2001a) (adjusted to the velocity offset of the O III A5592 hne that presumably forms 
deep in the photosphere and is less affected by expansion in the atmosphere) to find the 
peculiar component of space motion, Vp^c (following the methods described in Bcrger & Gies 
(2001)), and this quantity is listed in the final row of Table 2. We find that there is sat- 
isfying agreement between the predicted and observed space velocities, especially towards 
the higher mass range. Note that if the eccentricity had decreased since the the explosion, 
then the predicted velocities would be lower than the observed velocities, which is contrary 
to the results in Table 2. This strengthens our assumption that httle or no circularization 
has occurred in LS 5039 since the SN. 

The example of LS 5039 provides a strong confirmation of the predictions made about 
the eccentricity and runaway velocity that result from a SN explosion in a binary, and we find 
that both the observed eccentricity and peculiar space velocity can be consistently explained 
by our derived set of pre- and post-SN parameters. On the other hand, the system may 
be exceptional among MXRBs in the huge amount of mass lost in the SN. Indeed, LS 5039 
appears to be the fastest runaway object known among the MXRBs (Kaper 2001). The 
example of LS 5039 hints that there are other similar systems with low X-ray luminosity 
and small radial velocity variations, but this class of "X-ray quiet" SN descendants will be 
difficult to detect (Garmany et al. 1980). 

We thank the KPNO staff, and in particular Diane Harmer and Daryl Willmarth, for 
their assistance in making these observations with the KPNO Coudc Feed Telescope. Insti- 
tutional support has been provided from the GSU GoUege of Arts and Sciences and from the 
Research Program Enhancement fund of the Board of Regents of the University System of 
Georgia, administered through the GSU Office of the Vice President for Research. 
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Fig. 1. — The average spectrum of LS 5039 after shifting for orbital motion and convolution 
to the resolution of the 1999 spectra. The 2000 average spectrum is also plotted with the 
same offsets as the 1998 and 1999 averages {dotted lines) to emphasize the greater residual 
Ha A6563 emission observed in the first two runs. 



Fig. 2. — The predicted X-ray luminosity from wind accretion for an assumed 0-star mass of 
25Mq and a range of masses for the companion. The solid line illustrates the result for the 
average mass loss rate calculated for Mq = 25Mq (logM = —6.3 for mass loss in Mq y"^) 
while the dot-dashed hne and dashed line show the results for the strongest and weakest 
mass loss rates (from the 1999 and 2000 runs, respectively; logM is labelled in each case). 
The shaded region gives the estimated X-ray luminosity based on a distance of 2.6 kpc (the 
estimated distance for Mq — 25Mq) and the observed X-ray fluxes reported by Ribo et al. 
(1999). The best match occurs for Mx — 1.9Mq in this case. 



Fig. 3. — The mass plane diagram for LS 5039 with the constraints from the wind accretion 
model. The thick, solid line shows the estimated relationship between assumed 0-star mass 
and the mass of the X-ray source. The surrounding shaded regions show how the solutions 
change if we adopt the higher mass loss rate from the 1999 observations (lower Mx required 
to obtain the same observed X-ray luminosity) or the lower mass loss rate from the 2000 
spectra (upper Mx limit). Lines of constant orbital inclination are shown as thin, solid lines. 
The dashed line indicates the lower limit on Mx established by the lack of X-ray or optical 
eclipses (Ribo et al. 1999; Clark et al. 2001). 



Table 1. Range in Stellar Parameters 



Parameter Mq = IOMq Mq = 20Mq Mq = SOMq Mq = 40M( 



Ro {Rq) 5.2 7.4 9.1 10.5 

(km s^^) .... 2018 2400 2656 2854 

logM {MqJ-^) . -6.7 -6.4 -6.2 -6.1 

Mx (Mq) 1.2 1.7 2.1 2.4 

i?o(Roche) {Rq) 8.0 10.5 12.2 13.6 

d (kpc) 1.7 2.4 2.9 3.3 
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Table 2. Binary Parameters 



Parameter 


Mo = IOM0 


Mo = 2OM0 


Mo = 3OM0 


Mo = 4OM0 


jjinitial ^^"^ 


1.56 


1.56 


1.56 


1.56 


^initial 


14.2 


17.7 


20.2 


22.2 


^mitial 


5.7 


10.6 


15.2 


19.8 


AM2 (Mo) 


4.6 


8.9 


13.2 


17.4 


-Ro (-R©) 


5.2 


7.4 


9.1 


10.5 


i?o (Roche)""""' (Rq).... 


6.1 


7.7 


8.9 


9.8 


i?2(Rochc)""""' (i?o) .... 


4.7 


5.8 


6.5 


7.1 


Vsys (predicted) (km s~"^) 


120 ± 10 


154 ± 13 


178 ± 15 


197 ±17 


Vpec (observed) (km s~^) 


83 ±13 


118 ± 17 


145 ± 21 


168 ± 24 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

1998 
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